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Abstract 

The aim of this work is to compute the entanglement entropy of real and virtual particles by 
rewriting the generating functional of theory as a mean value between states and observables 
defined through the correlation functions. Then the von Neumann definition of entropy can be 
applied to these quantum states and in particular, for the partial traces taken over the internal 
or external degrees of freedom. This procedure can be done for each order in the perturbation 
expansion showing that the entanglement entropy for real and virtual particles behaves as In(mo). 
In particular, entanglement entropy is computed at first order for the correlation function of two 
external points showing that mutual information is identical to the external entropy and that 
conditional entropies are negative for all the domain of mo. In turn, from the definition of the 
quantum states, it is possible to obtain general relations between total traces between different 
quantum states of a cj/ theory. Finally, discussion about the possibility of taking partial traces over 
external degrees of freedom is considered, which implies the introduction of some observables that 
measure space-time points where interaction occurs. 


1 Introduction 

Entanglement entropy associated to a region of the space V has been extensively studied, where the 
degrees of freedom localized in that region are only taken into account and the rest is traced out. By 
using the von Neumann definition of entanglement entropy S = —Tr[/3ln(p)], it is possible to quantify 
the inaccessibility to the full system, p being the quantum state that results from the partial trace. This 
quantity has been widely used in several branches of physics, for example, quantum field theory (QFT) 
(m, 0,0, a and condensed matter physics and black hole thermodynamics (see m, m. 0 , 
m, BQ], m. |12) . |18) . |14] and [15) 1 and in particular for free quantum field theory with temperature 
(see [16], m and [H]), in curved space-time (see [19], [20], [21] and [22]), with excited states (see [28] . 
[21] and [2S]) and non-Lorentz covariant QFT [^. In the context of quantum field theory, geometric 
entropy of the free Klein-Gordon field has been related to the Bekenstein-Hawing black hole ([27] and 
m) and in general, free models can reveal features that are common to all quantum field theories, 
in particular, the interacting ones. In d dimensions it is shown that entropy behaves as a Laurent 
series starting in where e is a short-distance cutoff and the coefficients are functions on the 

boundary dV. The leading coefficient that multiplies to is proportional to the d — 1 power 

of the size of V, which is the area law for the entanglement entropy. In this work, the entanglement 
between real particles and virtual particles is investigated, that is, the subsystems considered are not 
partial traces over region of space, but over the intermediate states that are necessarily introduced in 
the perturbation expansion. In this work, the interacting (j)^ field theory will be considered and the 
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entanglement between external and internal propagators will be studied. Although the virtual states 
are a mathematical artifact of the perturbative expansion of the correlation functions, these states 
contribute to the physical mass, the vacuum energy and the coupling constant. On the other hand, if 
virtual states are not real or they do not exist (see [IH]); then it is feasible to trace out these states from 
the correlation function. In QFT, the particles that are created in these vertices are virtual particles 
because they are off shell; that is, they do not obey the conservation laws. In this sense, the conceptual 
meaning of the partial trace of the internal degrees of freedom is to neglect the virtual nonphysical 
modes. This is consistent with the experiments of scattering because basically what is measured are 
the in and out states. In turn, the interpretation of the integration of the internal vertices is to sum 
over all points where this process can occur (see [3D], p. 94). From the point of view of this work, the 
integration over the internal vertices reflects the fact that the virtual degrees of freedom are eliminated. 
Before computing the entanglement entropy, the quantum field theory formalism for a self-interacting 
system must be rewritten in a suitable way, which has been done in |31] and [32j with an application to 
nonrenormalizable theories in |33j . In m and [32] . the model has been applied to the renormalization 
of theory, showing that the renormalization procedure is equivalent to a projector that neglects the 
diagonal part of the quantum state defined through the correlation function, which is in turn a specific 
representation of the operator K. dehned in eq.(9.76) of [34]. In [33], the model has been applied to non¬ 
renormalizable theories, showing that renormalization group equation can be obtained. For the sake of 
simplicity, a short introduction to the main idea of papers |31j and [32j will be given in this section. In 
QFT, some (symmetric) n-point functions r(")(a:i,..., x^,) (like Feynman or Euclidean functions) can 
be considered; then the corresponding generating functional ([35], eq. (II.2.21), [3D], eq. (3.2.11)) can 
be defined as 


^ •■•,a;„)J(xi)...J(x„)]^(f'‘xi 

n=0^' i=l 


( 1 ) 


where 

t(”^(xi, = (fi \'p{xi)...(j){Xn)\ (2) 

and J{xi) are external sources. A convenient way to eliminate trivial contributions of single-particle 
propagators is by introducing a modified generating functional Z[J] for irreducible Green’s functions 
that is defined as 

W [J] = (3) 

The new generating functional Z[J] satisfies the normalization condition Z[0] =0 and it reads 


OO • ^ n 7~i 

iZ[J]='Y^'^— / T^'^\xi,...,Xn)J{xi)..J{Xn)\\d!^Xi 

n=0^' i=l 


(4) 


where in this case Tc^'^ {xi ,..., x„) are connected n-point functions that can be obtained by differentiation 


ri")(xi,...,x„) = ^ 




1 SJ{xi)...6J{Xn) 


|J=0 


(5) 


In turn, the connected n-point functions can be written in terms of the Lagrangian interaction density 
as (see eq.(II.2.33) of [35] ’H 

PC ^ 

Ti”^(xi, ...,x„)(^’^ = ^ / {no\T(l)o{xi)...(l)o(xn)£.°{yi)---C°{yp)\no)Y[d^yi (6) 

•' i=l 

Introducing ((DJ) in we have 


n -p 


= / {^o\T4>oixi):.(l)oixn)C°{yi)...C°{yp)\no) J{xi)...J{xn)Y[d‘^XiY[d‘^yi (7) 

n=0p=0^'7’’ i=l 


^In eq.© we have introduced the perturbative expansion of the correlation function, where the pi are the internal 
vertices. 
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This equation can be written as a mean value of an observable defined through the J{xn) sources in 
a quantum state defined by the correlation function (flo \T(j){xi)...4>{xn)C^{yi)---C^iiyp)\ flo)0 This 
procedure can be done for each correlation function of n external points. To define the quantum state 
we can consider some operator function F that depends on a set of vertices yi,...,yp and some new 
coodinates wi, ...,Wp in such a way that 


Y{yi,..,yp,wi,...,Wp)Y\5{yi - w,)= C°{yi)...C°{yp) 


( 8 ) 


where £^{yp) is the Lagrangian that appears in eq.Q- In [31] we have studied the theory for two 
external points and the corresponding operator can be represented by two different functional forms 


Fi{yi,..,yp,wi,...,wp) = Y[^(t)^iyi)(l){wi) 


P 


4! 


F2iyi,-,yp,wi,...,wp) = 




In both cases, eq.® holds. Then, inserting eq.® in eq. © we obtain 


[j] = EE 


n\ p\ 


(flo \T(j){xi)...(l){xn)F{yi, ..,yp,wi, Wp)| flo) 


n—Op—0 

J{xi)...J{xn)Y[Hyi - Wi)Y[d‘^XiY[d‘^yid‘^Wi 


n p 


i—1 2=1 2=1 

Now we can define two quantum operators in the following way 


gin,p) ^ J (^^Q\T(l){xi)...(j){xn)F{yi,..,yp,wi,...,Wp)\no) 

n p 

\xi,-,x^,yi, ...,yp) (xn+i, wi, ...,Wp\Wd‘^x^'^d'^yid‘^ 


2=1 2=1 


r P 

= / J(xi)...J{Xn)Y[d{yi - Wj) \xi,...,x^,yi, ■.■,yp) {x^ + l, ...,Xn,Wi, 
2=1 

n p 

Y\_d*Xii^d‘^yid'^Wi 

then, eq. can be written as 




i=l i=l 


n -p 




n=0p=0 

The quantum operator of eq. (na has the following form 




where 


/ n 

J{xi)...J{Xn) |a;i,...,X^) {x^ + l, ...,Xn\ Wd'^x^ 

2=1 


^In some sense, these observables will be recorded at the particle detector (see Ea, page 6, below eq.(2.6)). 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 

(14) 

(15) 


3 



and 


r(p) 



p 

- Wi) \yi, Vp) (wi,Wpl Wd‘^yid'^w^ = 

P 

\yi, ■■■,yp) {yi, ■■■,yp\W<Pyi 

2=1 


(16) 


is an identity operator acting on the yi vertices that appear in the perturbation expansion. The Dirac 
delta that appears as the coefficient of the identity merator can be considered as a particular choice of 
observable that physically implies no measurements The subscript ext in ea. (jl4p refers to the external 
points Xi and the subscript int to the internal vertices yi. Then, the generating functional of eq. © can 
be written as the mean value of the quantum operator Oext on the reduced operator as 

Tr{g^r.,p)oin,p)) ^ TriQ^foill) (17) 


where 


eilf = = J {yi, ■■■, yp\\yi, ■■■, yp) HAi = (18) 

2=1 

/ (Do \T(f){xi)...(j){xn)CPj{yi)...C°{yp)\ Oq) j \xi, ■■■,x^) (x^+i, f[d‘^x^ = 

i=i J i=i 

/ n 

T^'^^Xi,...,Xn) |a;i, (x^ + l, .■.,Xn\ 


where the subscript int refers to the partial trace over the degrees of freedom that represents the internal 
vertices of the perturbation expansion. This reduced state contains the divergences of QFT due to the 
correlation function (xi,..., a:„) that appears as the coefficient of the external reduced operator. 
What basically this model does is to duplicate the number of internal vertices. The new vertices are 
linked to the old vertices in such a way that identification one to one gives the Feynman diagram 
again. In turn, the way in which the observable in eq. m is written is suitable for a generalization, 
where the Dirac delta distributions are replaced by some nicer well-behaved distributions. This could 
be interpreted as if the interaction has been smeared out (see eq.©), although this implies that a 
nonlocal interaction has been introducedO This way of introducing the perturbation expansion of the 
generating functional of QFT allows to discriminate the internal vertices from the external ones in 
a simple way through the observable of ea. dUl) . If the quantum operator is normalized, the 

generating functional can be interpreted as a mean value between an observable defined by eq. (1141) and 
a quantum density operator defined in eq. (ED- This quantum density operator represent the probability 
amplitude for a particle to propagate from xi to X 2 , that is a distribution, that belongs to a C* algebra 
and where all the machinery of the algebraic approach of QFT can be applied to it. Nevertheless is 
not the purpose of this work to discuss the algebraic and analytical properties of these states, although 
some issues will be considered. 

As was shown in eq. CH), the mean value can be computed through a reduced quantum operator, 
where a partial trace over the internal degrees of freedom has been taken. This suggests that the 

^This point deserves more attention because a generalization of the Dirac delta can be introduced in such a way 
as to parametrize the partial trace, for example by using the Dirac delta representation ^(fc) = lim o or any other 

representation. 

'^An interesting point is that if the second operator function F 2 is considered, then the interaction term can be written 
as a mass term in the Lagrangian, where the coefficient is Ao^^(cd), which acts in another space-time point. This means 
that the mass is not constant and depends on what happens with the quantum field in another point where the interaction 
occurs. This point will be considered in future works. 
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partial trace over the external degrees of freedom can be considered as well. Both partial traces can 
be used to compute the von Neumann entropy defined as Sextjint = —Tr \_Qext/int^"^{Qext/int)\i where 
Qext/int ^re partial traces with respect the internal/external vertices respectively. This result could 
be important to understand the physical effects of virtual particles in real particles in the successive 
orders in the perturbation expansion and to obtain theoretical values of entanglement between these 
states. Although the procedure introduced in this section is simple, an important point has to be 
considered and is the ambiguity in the choice of the operator function F for a theory (see eq. Q ). As 
an example, for the case of the second order in the perturbation expansion with two external vertices 
and one internal vertex, the two operator functions imply the following open-loop connected Feynman 
diagrams 


Ql{xi,X2,yi,Wi) - 

^ (Ho \T4>{xi)4>{x2)Fx{yi,wi)\ Hq) " 

A(a;i - yi)A{x 2 - yi)A{yi - wi) 

(19) 

g2{xi,X2,yi,wi) - 

(Ho \T(j){xi)4>{x2)F2{yi,'Wi)\ Ho) - 

A{xi - yi)A{yi - wi)A{wi - X 2 ) 

(20) 


which represent two different Feynman diagrams, the first one with two external points and the second 
one with three external points. If we take the partial trace over the internal degrees of freedom we 
obtain the same reduced quantum state 

0exti^i^X2) Trintgi J d^yi {flo\T(l){xi)cj){x2)Fx{yi,yi)\no) A{ 0 ) j d'^yiA{xi - yi)A{x2 - yi) 

( 21 ) 

eextixi,X2) Trintg2 ^ J d^?/i {fio \T(l){xi)cj){x2)F2{yi, yi)\ flo) A{ 0 ) J d‘^yiA{xi - yi)A{x2 - yi) 

( 22 ) 

which is the first order contribution to the correlation function of two external points. This does not 
occur with the partial trace over the external degrees of freedom where two different results are obtained 


£’int(j/i:'*^i) ^ Txextgi ^ J d*xi (flo \T(l){xi)(l){xi)Fi{yi,wi)\ flo) ^ J d‘^xiA'^{xi - yi)A{yi - wi) 

(23) 

and 

glntivi^'^i) ^ Trextg2 ^ J d'^Xi {flo \T(j){xi)(j){x2)F2{yi,wi)\ flo) ~ J d^xiA{xi-yi)A{yi-wi)A{wi-xi) 

(24) 

In turn, the trace of both quantum states gives the same result (see figure [T|) 

Trg = Trext{Trintg) = J d'^xid'^yi {^o\T(t){xi)4>{xi)Fi/2{yi,yi)\^o) ^ A{ 0 ) J d*Xid'^yiA'^{xi - yi) 

(25) 

With these results, the only ambiguity is located in the partial trace over the external degrees of 
freedom which define the quantum state that represent only the internal vertices, which physically 
represent virtual propagation states. This could be related to the fact that in principle, the virtual 
states are an artifact of the perturbation expansion of the correlation function or perhaps that there are 
different ways of rearranging the internal vertices and links in such a way that the loop expansion of 
theory is obtained when the identification of internal and external vertices is done. This is particularly 
important when observables that depend on the quantum states are taken into account, for example 
the entanglement entropy. In this case, different values will be obtained for the external/internal 
entanglement entropies. In this work Fi will be considered only and entanglement and relative entropies 
will be computed at linear order in Aq. 

A different point of view for the idea behind the manuscript is that the correlation function of (j)^ 
theory is in fact the coefficient of a quantum density operator which is in turn a partial trace over 
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Figure 1: Scheme of partial traces over the possible quantum states. 


a larger quantum operator. This quantum density operator is the correlation function for a nonlocal 
interaction (j)^{yi)(j){wi) or (jP{yi)4P‘{wi). In the first case, the interaction is (f)^ which contains its 
own Feynman propagators. In other words, the following process can be considered for the operator 
function Fi: a particle in a definite momentum is prepared in the infinite past. When the interaction is 
turned on, this particle annihilates and two more particles are created. If the propagation of one of the 
resulting particles is not taken into account, then the propagator becomes a loop and the first order in 
the perturbation expansion for the correlation function of two external points is obtained. To ignore 
one of the particles is identical to putting two observables, one in the infinite past and one in the infinite 
future that measure plane waves, that has some definite values of the momentum operator, or in terms 
of group theory, two possible values of the eigenvalues of the mass operator in the reference frame at 
rest of the Poincare group, which is valid in the in and out states because the interaction goes to zero 
in those stages. One of the observables is the one which prepares the initial state and the second is the 
one that measures one of the particles that is a product of the interaction. In this sense, to take partial 
traces over the internal degrees of freedom is identical to make no measurement over the remaining 
particle, although it is there propagating in space-time. This point cannot clarify the discussion about 
the reality of virtual states, but it gives a physical reason of divergences in QFT: when intermediate 
states are not measured they must be traced out, then loops appear and infinities proliferate. In turn, 
the model introduced above is suitable for a generalization of observables, in particular, those that 
prepare and measure in and out states and simultaneously have an effect over intermediate states in 
such a way to avoid divergences. 

On the other hand, the partial trace over the external degrees of freedom can be interpreted as the 
following process: a particle is prepared in a definite momentum in the infinite past and in a particular 
space-time point, it annihilates and two other particles are created. If this product of particles is not 
measured but an observable that measures the exact space-time point where the interaction occurs 
is introduced, that is, where the initial particle annihilates in two other particles, then the resulting 
quantum state is that of eq. (1^51) . In the other case, where the trace is taken over the internal degrees of 
freedom, the observables only prepare and measure the in and out states, but in this case, one observable 
prepares a quantum state and the other observable measures the annihilation of the particle, which 
implies a measurement at one space-time point which would necessitate infinite energy. 

Another important point is to note that by applying a Wick rotation in the time coordinate, the 
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generating functional of correlation functions is identical to the partition function of statistical mechan¬ 
ics. Then it would be possible to study the entanglement of intermediate states, for tuq^ the Compton 
wavelength of the quanta which is the correlation length of statistical fluctuations and relate the results 
with those found in |38] , where area law for entanglement is obtained by considering partial traces over 
regions of space dictated by the decay of the correlations. 


2 (j)^ theory 


Before computing the quantum entropy of the reduced quantum states, we must take into account the 
algebraic structure of the Hilbert space. The quantum states can be written as 


= 


rr(p(")) L' 


,(».0) ^ ^(n,l) 



1 

Tr(p(")) 


+ CXD 

© g 

3=0 


(n,i) 


(26) 


where the supperscript n indicates the number of external points and i indicates the order in the 
perturbation expansion. The coefficient of each quantum state will be of the form 


...,a;„, 2/1, Pp, wi,Wp) = 

(Ho |T(()o(a;i)...())o(a:n)F(?/i,.., yp, wi ,..., Wp)\ Hq) 


(27) 


The trace readtH 

+ 00 

Tr(p(")) = ^(-*Ao)^'fH(„,,)rr(p(”’^)) (28) 

3=0 

where IT/n.i) is the weight factor (see |34) chapter 3) corresponding to the connected Feynman diagram 
and pi"’-!) is an operator that depends on the propagator of the respective Feynman diagram]^ The 
total quantum entropy can be computed as 


SM = _Tr[p(’^) ln(p(”))] 


(29) 


where S will be a function of Aq and some factor which will depend on the regularization scheme chosen. 
Up to first order in Aq, the quantum entropy in terms of p reads 


AobF(„p) 

(^(-.0))2fU(„,0) 


= ln(/3("’°)) - ln(p("’°))] 


0(Ag) 


(30) 


where = Tr(p*'"’*^). 


2.1 n = 2, zeroth order in the perturbation expansion 


In the case of two external points, at zero order in Aq, = rr[pg^’°^] and ln(p^^’°^)] must 

be computed and where VF( 2 _o) = 1 and W(^ 2 ,i) = 1/2- The quantum state at zero order is the free 
propagator 


( 2 , 0 ) 

Pext 


dip 

(27r)4 p2 _ ^2 


\xi) {X 2 \ d'^Xid‘^X 2 


(31) 


®It should be clear that the quantum states gl"! that depend only on the two external points are the partial traces 
over the internal degrees of freedom. 

®What this means is that the weight factor of the Feynman diagram which corresponds to the true perturbation 
expansion must be taken into account. For example, for n = 2 and 2 = 1, the quantum state is proportional to 
2^A(3 ::i — yi)A{yi — X 2 )A{yi — wi) which has a weight factor of 1/4 and the trace over the internal degrees of freedom 
gives a quantum state with a weight factor of 1/2, which is the corresponding tadpole diagram. The factor 1/2 is the one 
that must be take into account in eq. OH}. 
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Taking the Fourier transform by writing |a;i) = J e Igi) and (a; 2 | = / ( 92 ! the 

( 2 ) 

quantum state pg in momentum space is diagonal and reads 



(32) 

and the trace reads = i2TUAg, where 


' J (2^)4 (p2 - m2)^+i 

(33) 

and 2TV = f d^x = S‘^(p = 0) (see [30], page 96). Because is diagonal 

ln[p(^’°)] reads 

in the momentum basis. 

ln[pL’t^] - / /o u 1'4( 2 2 ) b) (P\ 

J (27r)4 p^ - mg 

(34) 

and Hp^ext^)] reads 


Trip^^ext^ ln(pi^’°^)] = -2TU (^Ag -P ixo) 

(35) 

where xo reads 

f d*p ln(p^ — mg) 

J (27r)4 (p2 _ ^2), + ! 

(36) 

Taking into account all the terms and using eg. (1301) at zero order 


Si^j, = ln(2rUAg) + 

(37) 


At this point it is crucial to compute Ag and xo with some regularization. Using dimensional regular¬ 
ization (see Appendix, ea. dMl) and ea. dM]) ! the entropy at order 0(Aq) reads 

sill = -^-l + ln(^^) + 0(e) (38) 

where e = d — 4 can be considered as a microscopic cutoff. In turn, the appearance of the logarithm of 
the microscopic cutoff e has been obtained in several works m, SO], ED, m and PH . The entropy is 
proportional to the dimensionless coefficient and reflects the fact that higher values of space-time 

volume or higher values of the mass of the propagating state increase the entropy. 


2.2 n = 2 , first order in the perturbation expansion 

In this case the total quantum state at first order in Aq and for n = 2 reads 


= J A{xi-yi)A{x2-yi)A{yi-wi)\xi,yi){x2,wi\d*yid‘^wid‘^xid*X2 (39) 

To compute the external entropy S'gj.) we can take the Fourier transform of eq. dMl) 

d^p 1 


P^exl^ = -*Ao 


\P) {P\ 


(40) 


J (27r)4 (p2 _ 

This quantum state is diagonal in the momentum basis, so ln(pg^’/^) can be computed as in the last 
section. Using eq. (sni) at order Aq, the entropy contribution reads 

XoAi' 


Aq 

2 


^ext ^ 2 

1 I 
-f 


dTT^e 


Xi 


27 — 1 -I- ln( 


TOn 


lOTT^/i' 


( 41 ) 


















where we have used / 3 *-^d) = —i 2 TyAoAi and we have introduced a mass factor to maintain the 
coupling constant dimensionless. From the last equation, the contribution at first order in Aq contains 
a microscopic divergence. Considering the last result and eg. (1551) the total contribution up to order Aq 
reads 


0 ( 2 ) _ 


1 J_ 

2 47r2 


. 1 , , . 1 , .TTIqTV Aq 

-1 -ln(e)-- + ln(^L^) + 


dTT^ 


327r2 


27 — 1 + ln( 


TOn 


167r2^' 


(42) 


This result implies that the reduced state at first order in Aq increases the entropy when virtual states 
are traced out. This point can be detailed as follows: if the quantum state of eg. (1391) normalized by 
dividing it by = —i2TVAoAi is considered only, then both quantum internal and external 

entropies read 


slip = -Tt[P1 HpH)] = ^ + In {2TVA,) = 


(43) 


4 

€ 


M^)+Oie) 


and 


s£" = -rrhg;‘>ln(p£‘>)l = f 


where we have used 


( 2 . 1 ) 


Pint 


= Tr,_ 


= - 1 + ln( 

e 


.( 2 . 1 ) 


47r2e 


f ln( 2 TyAo) = 
0{e) 


(44) 


d'^p 


^Tr(p(2d))J 2TVAo J {2tt)* - ml 


\P) {P\ 


(45) 


and we have used the results of Appendix A. The result of eq. (HI) is identical to the result of eq. dSIl) 
( 21 )... 

for the entropy of , which is a particular case of a general structure that it will be described in 
Section III. In turn, by applying the Fourier transform to the quantum state p^^^) t)e shown that 

it is not diagonal in the momentum basis 


.( 2 . 1 ) _ 


d'^pi d'^p2 d‘^p3 


2rCAoAi J (27r)'^ (27r)'^ Pi ~ 11*0 P 2 ~ P 3 ~ 


\Pl,P2+P3-Pl) {P2,P3\ (46) 


That is, the last equation implies a momentum entanglement as studied in |43j . Then, ln(p^^4)) cannot 
be applied unless the quantum state is diagonalized, then the von Neumann entropy must be computed 
in a different way. In particular, a family of functions called the Renyi entropies Sn where the limit 
n = 1 reproduce the von Neumann entropy is defined as 


S =-^HTr[p-]) 

and can be used to compute To compute the n—th power or p( 2 ,i) ji; be noted that 

I f d^pi d^P2 d‘^P3 1 1 V{P2+P3} 


(47) 


(p( 2 . 1))2 ^ 

where 

In turn 


(2TI4AoAi)2 J (27r)"‘ (27r)'^ (27i)^ Pi — Wq P 2 ~ ^0 Ps ~ ^0 

d'^pi 1 1 


i?(P 2 +P 3 ) = J 


(27r)4 (p2 _ ^2)2 +P3 -Pi)2 - mg 

d'^pi d'^p2 d‘^p3 1 I rf{p2+P3) 


\Pl,P2+P3-Pl) (P2,P3| 

(48) 

(49) 


{2TVAqAiY J (2t:Y {2t:Y {2ttY p\ — p\ — m^ pl — m, 


\pi,P 2 +P 3 -P 1 ) (P2,P3| 


( 50 ) 
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can be written as 


In this way, the n—th power or 


ip'wy = 


d'^pi d^p2 d'^Ps 1 


1 ??” {P 2 +P 3 ) 


{2TVAqAi )^ J (27r)4 (27r)‘* (27r)4 — TOq P 2 ~ P 3 ~ 

Taking the trace 


|Pl,P2 +P3 -Pi) (P2,P3\ 

(51) 




2TV 


_ f d*pi d‘^p 2 1 ??" (P 1 +P 2 ) 

{2TVAqAiY J {2ttY {2ttY {jpI - mlY pl-ml 


(52) 


Introducing the following change of variable pi + P 2 = f and d'^p 2 = d'^r, the last equation can be 
written as 


rr[(p( 2 ,i))"] = 


2TV 


{2TVAqAiY J (27r)4 


rf{r) 


Computing |jln(Tr[(p(2.i))n]) 


d , , f (2 u ln[77(r)]77(r) 

- HTriipi^^^^n) = - H2TVA,A,) + Y^r(r) 


(53) 


(54) 


where a derivative under the integral sign has been taken. Finally, taking the limit n —>■ 1 


where 


and 


Using eq.(8.19) of [34], r/(r) reads 


^(2,1) = ^+ln(2TUAoAi) = 


A = 


d\ 

(27r)4 


ln[? 7 (r)]? 7 (r) 


B = 


d'^r 

(27r)- 


-r/(r) 


(55) 

(56) 

(57) 


.r(3-i) 

7]{r) = —% 


(47r)'i/2 


/ da;(l — a;)[r^a;(l — a:) + tuq] 2 “^ = 
Jo 


(58) 


arctan( 


= —%■ 




%'K^r\J—rn^ — 


0{d - 4) 


and the A and B coefficients have been computed in Appendix B.Taking into account the results of 
eg. (1551) . ea. dTSl) and ea. lHTl) . mutual information can be computed and reads 


/(p(^’i)) = + ^(p 5')) - 5(p(24)) = ^ _ I + in(2TU) (59) 

=-hi — r + ln(2TOQTU) 

In figure m the finite part of total entropy entanglement entropy and \ mutual 

information /(p^^’^l) and the sum of external and internal entropy is plotted as a function 

of bare mass mo using that TV = 1. The subadditivity of a bipartite system (see HD) can be seen 
to be obeyed by noting that the dashed curve corresponding to is larger than 

In turn, entanglement entropy and total entropy go as ~ 41n(mo) as occurs for entanglement entropy 
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Figure 2: Total and entanglement entropies, mutual information for the first order perturbation expan¬ 
sion of the two correlation function. 


for different geometries (see |45jl. Mutual information and external entropy are similar only for low 
values of mg, where a difference can be seen in the inset of figure [5] Internal and total entropy behaves 
similarly for larger values of mg, although differences can be obtained for values of mg near mg = 1. 
These results imply that the partial trace over the internal degrees of freedom at first order in the 
perturbation expansion has a negligible effect on mutual information and the total entropy of the 
system. Of course, as is expected, entanglement entropy for external and internal states are larger than 

total entropy, although it is bigger for than . In turn, the finite part of the conditional 

(2 1 ) (2 1 ) 

entropy for and can be computed and both results do not depend on mass mg 

- S{ptl^) = r + 1 - ln(87r2) ^ -0.102 (60) 

and 

= r - 2 - ln(87r2) ^ -3.102 (61) 

Both results have negative values for the whole domain of mg at first order in the perturbation expansion 
which reflect the fact of the quantum nonseparability of the total system and the coherent information. 


2.2.1 Nonperturbative approach 


General properties can be obtained by not taking into account the perturbative expansion in Ag, but 
considering the spectral representation of the two-point correlation function (17 |T^(a;i) 0 (x 2 )| 17) where 
17 and (/)(a:) are the vacuum state and the field operator of the interacting theory. Following the 
procedure of the Introduction, a quantum state, which is the partial trace over the internal degrees of 
freedom can be defined as 

= - / (17 |T(()(a;i)(^(a;2)| 17) Icci) (a;2| (62) 

V J 


where rj = reads 

ry = y (17 |T(^(xi)(()(xi)| 17) 


(63) 


and is introduced in the definition of the quantum state of eq. (Ea to have a normalized quantum state 
Tr{Q^^'>) = 1. Using eq.(7.6) of [30] we can write the last equation as 



+ 00 


27r + ie 


|xi) (X 2 I d‘^xid^X 2 


(64) 
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where a{M^) is a positive spectral density function (see eq.(7.7) of [30]) which contains one-particle 
and multiparticle states. Taking the Fourier transform, the quantum state in momentum space is 
diagonal and reads 


Q 



d'^p f+°° dM^ , , , , 

(2^F i + <*’1 


(65) 


then 


where 


and where 


7 = 


r] = i2TV"f 
'■+°° dM^ 


2tt 


a(M")Ao(M") 


Ao(M 2 ) = J ■ 


d^p 


1 


(27 r )‘* _ ]\/[2 _|_ jg 

The entropy of the quantum state defined above reads 


S'(2) = ln(i2Ty7) - - 


d'^p 


^+oo 


dM"^ ia{M^) 


■In 


f*+oo 


ia{M^) 


27r p^ — - 1 - ie iJq 27r p^ — - 1 - ie 

Considering only the contribution to the spectral density of the one-particle states 

cr(M^) = 2ttZ 6{M'^ — m^) + multiparticle states 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


where Z = |(fl |^( 0 )| Ao)|^ is the field-strength renormalization and m is the physical mass, the entropy 
reads 

5 = ln[2TyAo(m2)] + |^ (71) 

where Ao(to^) and Xo(w^) are the functions defined in ea. (l551) and ea. (l551) with rug replaced by the 
physical renormalized mass m?. This last result is in fact identical to the result obtained in eg. (1571) . 
The first integral of the last equation is the known free propagator and the second integral has been 
computed in Appendix A. Although eq. (ED does not contains any approximations, its depends on m 
which is a renormalized constant that depends on Aq due to the renormalization group equation. 

Finally, as was said in first section, by considering O G where B is the set of all bounded 

operators that form an algebra A acting in a Hilbert space "H, a linear form p over this set can be 
defined. In particular, the Gelfand-Naimark-Segal theorem m implies that each positive linear form 
ip has a representation in the set of bounded operators A. This allows one to define a scalar product 
on A as ip{0) = Tr{'K{ip)0) where 7 r((/?) = p is the representation of ip over A. For two external points, 
the quantum state can be used and the trace with an observable defined as 


O = 


J{xi)J*{x2) |a;i) {x2\ d‘^xid‘^X2 


(72) 


reads 


Tr{p^^'>0) = 


1 


/*+oo 


dM^ 


p{M^ 


[ d*p 

Jip) 

2 

j (27r)4 

- -I- ie 


(73) 


2TVj Jo 2 tt 

where J{p) is the Fourier transform of Jix) and 7 is defined in ea. (l67L which is the normalization 

2 

factor introduced as having Tr(p(^)) = 1. In the case of plane waves, J{p) =1 and Tr{p^‘^^0) = 

{2TV)~^ > 0. It is not difficult to show that ip{0*0) = Tr{pO*0) > 0 and that <p(0*) = Tr{pO*) = 
Tr{pO) = ip{0) where the top bar indicates conjugation (see page 122 of [35]). In this sense, the 
positive linear form can be considered a state because <p(/) = Tr[p) = 1 and all the machinery for 
entanglement entropy can be applied. 
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2.3 n = 0, first order in the perturbation expansion 

In a similar way, we can compute the total entropy of the quantum state related to the vacuum-vacuum 
amplitude. We can consider the inner product between (n| and |n) (see page 87 of [30] or eq.(29) of 


m) 


where 


(fl I n) = 1 = 


e*‘®“^^Tr(p(o)) 


Tr(p(°)) = ( flo 


exp[—z f dtHj{t)] 
J-T 


—zAc 


-T 

2 



Sin ) — 1 


(fyi (Slo |(('o(2/i)| ^lo) 


(74) 


(75) 


d'^yid!^y 2 (flo |<('o( 2 /i)')>o(j/ 2 )| flo) 


where we are considering that the quantum states up to first order in the perturbation expansion reacfl 

1 


pi.m = 


2TV 


\yi) {wi\(ryi(rwi 


(76) 


1 


= ^ / d'^yid'^wi (Slo |<)>o(2/i)^o('u;i)| Slo) |yi) (wi| 


In this case, the quantum state, which represents the generating functional of the n = 0 external points, 
contains no real particles; then, the entropy will be related to the process of creation and annihilation 
of virtual particles. Using eq. (PI and eq. dZSI, the contribution at Hrst order in Aq to the entropy reads 


= ln(2rU) 


1 - ^TUAo 


An 


which, by using dimensional regularization, becomes 

^(0) 


= 1 _ ^TUx 
ln(2rU) 4 ^ ^ 


(77) 


(78) 


( < -2 

■ 

\647r4 

647r^ 


7-l-ln(^)) -87r2 


-1 


1 


I5367r4 


[d -t Bm'^ + + 0(e) 


where 

and 

and 


A = 967r^[l — 7 -I- ln(47r^^)] 

B = 18 -|- 12(7 — 2)7 + Tt^ + V 2 ln(47r/i) [ln(47r/i) -|- 2 — 27] 
/(toq) = 48 TOoln(mo) 


-1 


.In(^) 

^ Any. 


(79) 

(80) 

(81) 


In hgure[3l the finite contribution at first order in the perturbation expansion for the entropy of virtual 
processes is plotted against toq for different values of mass factor /i, which is a renormalization scale, 
that is, the renormalized couplings depend on the y value in such a way as to obey the renormalization 
group equations. Because ztzq depends on then the entropy computed in eg. (1331) and eg. (1751) depends 
on 2TV and the mass scale y considered. In all the cases, entropy has a minimum for a particular 
value of ztzq which decreases when y increases. It should be clear that the lack of knowledge that 
leads to the entropy contribution at first order for the zero point correlation function comes from the 
absence of measurement of the propagation of a real particle from one space-time to another as can 


’^Eq. II74II implies that Tr{p^^^) = KHo I 
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Figure 3: Entropy of the zero point correlation function at first order in the perturbation expansion for 
different values of mass factor /i. 

be seen from eq. (EH). This measurement is not done at all because what is being considered is the 
process of vacuum to vacuum amplitude that occurs in short periods in time. In this case, the quantum 
entropy computed is not important. But it is useful in the case in which the virtual process converts 
to a real process, for example see m, where the effect on the entanglement entropy of real particles 
coming from a common virtual pair is considered or the opposite process where a massive pseudoscalar 
particle decays into a particle-antiparticle pair (see |48j 1. For example, if the first order contribution to 
the vacuum correlation function is considered and one of the loop propagators is cut, then we obtain 
the partial trace over the internal degrees of freedom of the first order contribution of the two-point 
correlation function. If it is possible that some physical process that involves these two quantum states, 
for example a loop propagator converting in a real particle propagating in space-time then the entropy 
of both states can be computed separately as it was done in Section II and then compared. 

3 General considerations 

A relation between traces of and p^^'> can be found by noting that the trace over external 

points in or gives contributions to the perturbative expansion of As an example, the first 
( 2 ) 

order contribution of p\J^ can be considered 



(82) 


and we take the trace on the reduced state we obtain the trace 


Tr{p^^’^'>) = TrextipfxP) = -*AoAo / A(a;i - yi)A{xi - yi)d'^yid‘^xi 


/ 


(83) 


which is in fact one of the second order contributions to Tr{p^^^) divided by Aq (see figure IH blue 
boxes) 



that is 



(85) 


In turn, the other contribution to Tr(p^°’^)) 



(86) 
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Figure 4: Scheme of partial traces over the possible quantum states. 


is in fact proportional to the total trace of (see figured black boxes) 


= -zAo J A‘^{xi - yi)(fyid^xi 


(87) 


that i^ 


Tr(p(4.i)) 


—iAo 


( 88 ) 


In a similar way, a relation between the connected Feynman diagrams of n = 2, n = 4 and n = 0 
correlation functions can be done, showing that a general expression can be obtained for a (j)^ interaction 
in terms of the quantum state traces 


rr(p(°)) = 1 - zAo[Tr(p(4)) + AoTrip^^^) - Ag] (89) 

where the product of loop propagators Aq corresponds to the trace of a free propagator multiplied by 
Aq, which is the zero order in Aq of the n = 2 correlation function (see figured the first Feynman 
diagram of Tr{p^'^'>)). The last equation can be related to the generation of correlation functions from 
vacuum diagrams (see Section 5.5 of [31], page 68), where for example, by cutting one line to the first 
order vacuum diagram we obtain the first order contribution to the two-point function. In terms of 
quantum states, the ’’vacuum to two-point or four-point functions” way cannot be done because cutting 
a line implies introducing a new propagator, which implies the introduction of at least two quantum 
fields that act as external points in the correlation function and this, in algebraic terms, implies to 
introduce a new Hilbert space in the algebraic structure. Then, because the enlarged quantum state is 
not a tensor product of quantum states of each Hilbert space, then there is no operation that allows 
finding the correlation functions from vacuum diagrams. The opposite way is the one introduced in 
eg. (1331) to ea. (l88L where by using correlation functions it is possible to obtain the vacuum diagrams. 
In fact, because the procedure is the same and because the weight factors in each diagram match as 
it is pointed out in [34], then the weight factors match as in eas. (l82ll8'8l) . which enables us to obtain 

®It must be said that the trace over pHl is taken by considering identical arguments, that is Tr(p^) = 
f(z,z\p(^)\z,z)d^z. 
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eq. dHU). This equation can be generalized for a general (jf interaction: 


-—1 

Tr(pW) = 1 - iM-C- - 1)A| + ^ (90) 

j=o 

which puts all the <j)^ theories on an equal footing. From ea. dMl) . eq. (1711) and ea. (|55)) . Tr{p^'^^) can be 
obtained as 


Tr{p^^^) 


|(flo I 


1 + Ao(m2) [Ao(mg) - 2TVZAo{m^)] 


(91) 


where m§ in turn can be written in terms of bare mass In eq. (4.538) of |49) . the general relation 
between and mg is shown where the coefficients of the expansion depend on the renormalization 
prescription, for example the mass-independent prescription (page 137 of [H]). In this prescription, 
mg, Ao and Z depend on the renormalized coupling constant A, so then the r.h.s. of the last equation 
can be written as a function of A, the vacuum energy Eq and the space-time volume 2TV. 


4 Conclusions 

In this work, the entanglement entropy between real and virtual propagating states has been computed 
by rewriting the generating functional of the (/>^ theory in such a way to separate the internal propagators 
from the external ones by considering the initial and end points as labels of position states belonging to 
particular Hilbert spaces. This formalism does not introduce new mathematical elements, but allows 
us to introduce the idea of observables that measure the in and out states and do not measure the 
intermediate states. In this way, the divergences of quantum field theory appear from the partial traces 
taken over the internal vertices of the respective Feynman diagrams. Entanglement entropy has been 
computed for two external points at first order in the perturbation expansion. It is shown that mutual 
information between real and virtual particles increase with bare mass and that the conditional entropy 
is negative, which implies that virtual and real particles are highly entangled. In turn, the entropy of 
a virtual process such as a vacuum to vacuum amplitude has been computed showing that minimum 
values are found for particular choices of mass factor p. Finally, general results can be found for total 
traces of quantum states that represent zero, two and four external point correlation functions. These 
results imply that the formalism introduced in this work naturally relates different correlation functions 
for a general (jf theory. 
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A Appendix A 


In this appendix Xj is computed using dimensional regularization, where Xj is that of ea. dMl) . In the 
hrst case 


Xj = 


d^p ln(p^ — TOq) 
{2t:Y (p^ — TOQ)t+i 


(92) 


By applying Wick rotation to perform the four-dimensional integral in four-dimensional spherical coor¬ 
dinates, where po is switched with ipo, then = —p|; and by using ~ and introducing 
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the following change of coordinates x = the last equation can be written as 


Xj =■ 




d-2{i+l) i-l 

0 / 


(47r)'^/2r(d/2) 7o 


1 / ^6 \ 7 — - 

lll(- -)X^ 2 


(l-a;)2 ^dx 


(93) 


The result reads 


Xj = 




m, 


d-2{j+l) 


(l)^r(j + l-d/2) 


(47r)<^/2 


r(j +1) 


Hj - +\n{-ml) 


(94) 


where Hj is the Harmonic number of order j. In a similar way, the functions Aj can be computed for 
an arbitrary dimension d. In several textbooks these functions are shown (see Appendix A.4 of |30] 1 


^j — 


d‘^p 


r(j + 1 - d/ 2 ) 


{2 tt )'^ { p ^ — TO§)t+i (47r)'^/2 


r(j + 1 ) 


(95) 


These results will be used in the main sections of the manuscript. 


B Appendix B 


In ea. (l55ll . A and B coefficients as a function of mo must be obtained. To do so, the function T]{r) has 


been obtained in eg. (15^ . Then, by introducing the following change of variable, s = 
can be written as 

where 


iy/m'^+r 


B = bm-Q 


eq.dSZl) 


(96) 


b = 


1 f arctan(s) 

Stt^ Jq (1 + s2)2 


ds 


(97) 


where the last term is is not bounded for x ^ —i. In a similar way, by using the same change of 
variables, the A coefficient reads 


where 


In this way 


I 


a = 


A = mg[a + 51n(TOQ)] 


* arctan(s) A(1 + s^) arctan(s), 
ln[-i-y:- ]ds 


StT^ _/q (I + StT^S 

A mQ[ai + 51n(mg)] a 


B 


bm/i 


ln(mQ 


(98) 

(99) 

( 100 ) 


where a/b is not bounded in the limit s —>■ —i. Nevertheless, this divergence can be grouped with the 
microscopic cutoff divergences in the result of eq. (ESD. For the s —>■ 0 limit 


lim— = T = -[—27 — ln(4) + 12 + 3<C(3)] ^ 3.2663 
s->o 6 4 


( 101 ) 


where C(-2:) is the zeta function. 
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